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We analyze properties of degree and clustering of a hyperbolic geometric model of complex networks in 
small parameter case T < 1,20 < 1. We find that the probability of k-degree goes to 0 and the global 
clustering coefficient goes to 0 in probability too as the number of nodes N — œ for some specific growth 
R(N) of the region radius. Here the scale-free degree is failed and the connection between neighbors are 
very weak. The transition of properties of the model with the parameter o changes seems to show that the 
mobility is important to keep society full and stable communication, otherwise a silence society. Some 
analysis technique and method are first applied for such model. 
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1. Introduction 


Many models of complex networks have been proposed to capture some common properties of real 
world. Perhaps the most well known is the Barabasi—Albert model of preferential attachment [3], which 
is based on two generic mechanisms: (i) networks expand continuously by the addition of new ver- 
tices, and (ii) new vertices attach preferentially to sites that are already well connected. Their model 
reproduces the observed stationary scale-free distributions and indicates that the development of large 
networks is governed by robust self-organizing phenomena that go beyond the particulars of the indi- 
vidual systems. Another common property of social networks is the cliques form, representing circles 
of friends or acquaintances in which every member knows every other member. This inherent tendency 
to cluster is quantified by the clustering coefficient or the global clustering coefficient. Although the 
Barabasi—Albert model is effective in capturing the scale-free degree sequence, but fails to generate 
networks with the clustering phenomena. There are growing indications that the tendency for affiliation 
of neighbors of the node is a manifestation of the network having a geometric structure [11], based 
upon the notion of a geometric random graph. A statistical mechanical models of complex networks is 
introduced in [14], and, in particular, it is applied to random geometric graphs in hyperbolic space [12]. 
The model employs a random mapping of graph nodes to points in a ball centered at the origin in the 
hyperbolic space H¢*!, giving these nodes hyperbolic coordinates, and connecting probability about the 
hyperbolic geometry. Therefore the random geometric graphs in hyperbolic geometry provide extremely 
promising models for the structure of complex networks. 

In [16], we extend the hyperbolic geometric random graph model of Krioukov et al. [12] to arbitrary 
dimension and considered this model in five regions of the parameter space for any dimension, and 
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we found the power law (scale-free) expected degree distribution for two regions, and for others, the 
power law disappears and the probability of k-degree goes to zero. We remain T < 1 and 20 < 1 case 
no research since the main difficulty lies in that the probability of connection between a fixed node u 
and a random node v is mostly contributed by ry < R—r, + @(N) and the previous explicit formula of 
hyperbolic distance is powerless to it. 

In this article, we continue our study for the remaining parameter region T < 1 and 20 < 1. We 
explore the hyperbolic distance formula deeply and overcome some technique difficulty, and got the 
estimation for the probability integral of connection, hence we may analyze the degree and clustering 
property of the model. We find that the probability of k-degree goes to 0 and the global clustering 
coefficient goes to 0 in probability as the number of nodes N — œ for some specific growth R(N) of the 
region radius. Some analysis method and techniques are first applied for such model. For the clsutering 
study, we borrow partly some method of Fountoulakis for tT < 1,20 > 1 and d = 1 in [4]. 

In generalization of analysis results in [4, 5, 16] and this article, for N nodes random geometric 
graphs G(N,¢,t,0,V,d) in the hyperbolic ball with the graph Hamilton energy parameter T < 1 and 
the rescaling hyperbolic ball radius Zy = In y we have the classification as follows. There exists the 
power law of degrees for the 20 > 1,d > 1 in [5, 16]; and there don’t exist scale-free degree distributions 
and the probability of k-degrees go to zero as N — œ for the 20 < 1,d > 1 in [16] and this article; the 
global clustering coefficient will go to some nontrivial constant in probability for the 20 > 2,d = 1, 
which mean the connection ratio of neighbors are scale-free in probability, and the global clustering 
coefficient will go to zero in probability as N — œ for the 1 < 20 < 2,d = 1 (we believe that such 
results are right for d > 1) in [4]; we prove that the global clustering coefficient has the same to-zero 
tendency in probability for the T < 1,20 < 1,d > 1 in this article, which mean that the connection of 
neighbors will become less and less possible. We don’t study the global clustering coefficient for the 
20 = | case since plenty of computing, but we believe that it is the same to-zero tendency. Altogether, 
we can see that the density function (1.3) of 206 = 1 is exactly a critical probability for the scale-free 
degree in the random graph and so is o = | for the clustering as Albert and Barabdsi mentioned in [1]. 
There are also many papers such as [9, 10, 12] to study the hyperbolic geometric random graph model 
with numerical method, which approximating results prompt the development of the analysis method. 

Why the degree and clustering of networks become more weak for the growing speed of the radius 
Ay = In¥ as the parameter o become smaller in the hyperbolic geometric random graph model. We 
may observe that the probability of nodes (1.3) staying at the interval [r,r + dr] of the small radius r is 
increasing when the parameter o become smaller, and correspondingly nodes are more likely to stay at 
the closer region to origin, for example r < 5 with d = 1, the density of this region is adding when the 
parameter o become smaller, hence two nodes in this region possibly are closer and easier to connect 
since the formula (1.2); on the other hand, we can observe that the probability of nodes lying in the 
close region to origin is decreasing when the hyperbolic ball radius is increasing, for example r < a 
with d = 1, the probability of nodes staying at this region is decreasing as the radius growth. Hence the 
two opposite tendencies interact each other: the willingness of nodes staying at the interval [r,r + dr] 
of the small radius and hence dense as the parameter o become smaller is resisted by the radius growth 
at the speed Zy = In x, in contrast, the willingness of nodes staying at far region away from the origin 
as the parameter o become larger are consistent with the influence of the radius growth. Finally there 
achieves a critical point at 20 = 1 for the scale-free degree and at o = Ifor the clustering property of 
networks at the fixed growing speed of Zy = In® . This phenomenon in the model is more like the 
principle in human society, the willingness of nodes staying at the interval [r,r + dr] of the small radius 
are representing the willingness of people’s immobility and keeping at small circle, otherwise mobility. 


If people are dislike to move, the communication among people will become less and less when the 
society expands as some kind speed of population, and it finally form a silence society, on the contrary, 
if the likelihood of mobility of people is increasing, the communication among people will become 
more and more easy and frequent, so mobility makes the society full and stable communication, i.e., the 
scale-free will occur. 

As the above analysis, should the mobility be another generic mechanism for networks, especially 
social network? 


1.1 Model Introduction 


We analyzed a class of exponential random graph models in [16], in which a random graph G = 
G(N,C,7,0,V,d) with N vertexes, denoting the vertex set by Vy, has the positive model parameters 
Č, 0, v, and T, whose elements are randomly distributed into a ball of radius R centered at the origin in 
the hyperbolic space H@*!, integer d > 1, with a probability density 


pu(%) = p (r)po (6) = p (r)p1 (01) --- Pa( 8a) , (1.1) 


and the probability of occurrence of an edge between vertices u and v is 


1 
Puy = — nR s (1.2) 
l+e T 
Where x = (7,9) is the spherical coordinates on H?+! with the usual coordinates 6 = (01, 02,--- , 04) 


on S¢, O; € [0, m) for k =1...d—1 and 04 € (0,27), and r € [0,R). The density functions p(x) = 
sin?~* O /Ia, with Ia = J sin? “0 dO for k=1...d —1 and pa(04) = 34, and 


= sinh? (o¢r) 


p(r) F (1.3) 


with Cy = fr sinh? (o Çr)dr. The distance function d, is the hyperbolic distance between two points 


> 


u= (Fy, ð) and v = (ry, 9,), which is given by the hyperbolic law of cosines 
cosh Édu y = cosh Gry cosh Er, — sinh Cr, sinh r, cos Ou v, (1.4) 


where O„» is the angular distance, i.e., the relative angle, between 6, and 6, on S¢. The parameter 
o governs the radial node distribution, v the node density, T the ‘temperature’, and —¢? governs the 
curvature of the hyperbolic space respectively. The conditional probability that a node u on the fixed 
position of the hyperbolic ball connects to a angular random node v at fixed radius r, is 


T sind! @ 


Tq, 


T 
Puy =) Pi (Ou) Puy d Ou» =} Puyvd yy , (1.5) 


which we call the angular integral. We let J,,,, be an indicator random variable which is 1 when there 
is an edge between u and v in the graph, and 0, otherwise. So the conditional probability that a node u 
on the fixed position of the hyperbolic ball connects to a random node v is 


= R pt oind-1 ra) 
Priluy = Iru, (JA =) | ice Le dO, ydry. (1.6) 
o Jo Ta 


The integrals in (1.5),( 1.6) are both only dependent on the relative angle 6, ,,, so the probabilities 
are both independent of the specific angular position 0,. We will omit this specific angular position 
sometime. The expected degree of the node u on the fixed position is 


(ku) = (k(ru)) = (N — 1) Prily = Uru] - (1.7) 


If we average the random movement of the node u, i.e., computing the node expected degree, we get the 
mean degree 


W = [wp rar. (18) 


We also let D, = Yrevi Iu» is the number of connections to the vertex u, i.e., its degree, where Vy = 
Vy \ {u}. We called the connection number equaling to the positive integer k by k-degree, and its 
probability denoted by Pr[D,, = k]. Throughout we let R = R(N) > œ, @(N) =0(R(N)) > œ, o(1) +0 
as N — co and op(1) means — 0 in probability as N — œ. We also use the Big Theta © notation for 
some function f(x), i.e., O(f(x)), which means, Sc; ,c2 > 0, s.t. 


cif (x) < OF (x)) < cof (x). 


In this article, the function f(x) is related to N and N — œ, so we always assume Big Theta function 
©(f(x)) true at N enough large situation. To simplify the notation we will rescale variables as follows 


d d 
(.&u,1) = 2 (RT), oi (N) = Saw), 
There are also some useful computing results: Cy = (1 + 0(1)) a and 
r sinh? z 
f an dr = (1+0(1))e%-®) = (1 +.0(1))e2°(-##) (1.9) 
0 d 


20 -2n 
for r > 1, and p(r) = (1 to(logae r from [16]. There have some notations about the 


clustering in the network model: 


e A(u,v;w) represents the event that the triple u,v, w forms an incomplete triangle pivoted at w, in 
other words, u,v, w forms a path of length 2 with w being the middle vertex. 


e A(u,v,w) represents the event that the triple u,v, w forms a complete triangle. 


e T =T(G),A = A(G), the number of complete triangles or incomplete triangles in the graph G 
respectively. 


1.2 Main results 


We study the degree and clustering property of the random graph G = G(N,¢€,7,0,V,d) with N vertexes 
and parameter T < 1 and 20 < 1 in arbitrary dimensional hyperbolic space H¢*!. The main results of 
this article are as follows. 


THEOREM 1.1 Lett<1,20<1andd È 1, then 


E R 
Prifiy = 1 ru, 94] = [ PuvP(ry)dry = O(e-70™) (1.10) 
0 


for any r,. Mean degree 
(k) =@ ((N- 1) Bue?) (1.11) 


The probability of k-degree 


as N — oo, for any given k > 0, if let Zy = në. 


THEOREM 1.2 Lett < 1,20 < 1 and d È> 1, then 


B(a)=3( Jo (e274), (1.12) 


There exists the the lower bound estimation and the order estimation for E(T) 
N 302, 
E(T) 2 e H and E(T)=o(E(A)). (1.13) 


THEOREM 1.3 Lett < 1,20 < 1 and d > 1. Also let Zy = In. The global clustering coefficient 
C2(G) of the graph G(N,¢,7,0,V,d) will go to zero in probability, i.e., 


3T(G) _3E(T)(1+0p(1)) 
A) = Te = Faro) 


+0, 
as N —> œ. 


The paper is organized as follows. In Section 1.1 we introduce the model. In Section 2 we give a 
useful preliminary result. We give some degree estimation and prove Theorem 1.1 in Section 3. We 
discuss incomplete and complete triangles and prove Theorem 1.2 in Section 4. Finally we estimate the 
variance and prove Theorem 1.3 in Section 5 and 6. 


1.3 Simulations 


We adopt R software https://www.r-project.org/ for statistical simulation. According to the law of large 
numbers in probability theory, we use the mean to approximate the expectation in R program. We simu- 
late the connecting probability of a fixed node in two kinds of dimensions, k-degree and global clustering 
coefficient of incomplete and complete triangles. One also can use some usual ComplexNetworks 
toolkit within Cactus [2] for the simulation. 


2. Preliminary lemma 


The hyperbolic distance d, in (1.4) can be explicitly expressed according to the following lemma. 


LEMMA 2.1 Let h,hı : N—R® such that A(N), hı (N) — œ% as N — œ. Let u,v be two distinct points in 


: A F : A 1 : 
H¢+! with Ou» denoting their relative angle. Let also Ou v := (e726 e726) 2, If we assume u,v with 
Tu, ty > h(N), then as 


we have 


2 hg 6.) 
duv = Tu + ry + č Insin (% ) +0 ( anp) À (2.1) 


here |O (( funy? ) |<3a 3 few )? uniformly for all u,v satisfying the above condition. Note that we can 
take hı (N) =h(N), for example. 


Proof. For ry,ry > h(N), we have Ô. » <S V2e~$"\). Thus, we may choose hı (N) to make hy (N) Ôv < 
T, hence the condition h (N ) Oy, < Oy) < 7 is well defined. Obviously, we can take hı (N) =h(N) so 
that A(N) 6, < T. Next, we prove (2.1). The right-hand side of (1.4) 


eb (Mutt) sin? 6. (: N (1+cos Guess ) into) 


cosh Cr, cosh Cry — sinh Cr, sinh Cr, cos Oy) 


2 2 2 sin” Puy 
C(rutry) 0, 6 
e . 2 Guy UV \2 
1+ ; 2,2 
5 sin’ ( Og) ) (2.2) 


Here, 0 < O(( a )?) < 2n7( i )?, and we also make use of Young’s inequality to get 


on Slur) = 926M <b (alr 4 e727) < G2 


u,v * 


Thus, 


2 x 
1 u,v AV 
cosh €r, cosh €r, — sinh ¢r, sinh Gry cos, > (š ) (: | oq: ») ; 


so from (1.4), we get 


1 u,v 1 
eda > LEL > wl). 23) 


Furthermore, from (1.4) and (2.2), we have 


Cdyy +In(1 +6725) = Ẹ (ru +r) +Insin? Pur +In ( oÊ -) ») ; 


and, from (2.3), we have the result 


2 u,v Â, v 
duy = (ru +rv)+ i a LO (6 7 r) ; 


where |O (Ge) |< 3m? (gu). 


Especially, for all u,v with ry +r, — R > @(N) with @(N) = o(R(N)) > œ as N > œ, we may choose 


(2.4) 


if 


4 
m A v2 30N) 
then (2.1) is achieved, since r,,7, > @(N) and 6, > hı (N)ĝ, y, where hı (N) = 2 oy 


oy further we 
have 


. Buv buy 
dy.y—-R rutryR+ F Insin( 5" )+O(( ge )?) 
T 


e T = é 


by 2 


o(( ie 
rutry-R , d 6, v (( u,v ") 
T T T 


NQu+nv-2y 
t 


with 


sO 


dyy—R = Nu+Nv-ZH . d 
l+e`T 1+Cwe T sint 
3. Degree estimation 


In this section, we prove Theorem 1.1. 


3.1 The connecting probability with a fixed node for arbitrary dimension 


We compute Pr[J,y = 1|ru, 6,,] in this subsection. 


=> 


3.1.1 Prilyy = 1|ry, 8u] on the condition r, > @(N). As the definition 


7 R 
Priluy = Ary, 6,,] = [ PuyvP (r,)dr, (3.1) 
R—r,+(N) R 
= I Pu,vP (r,)dry + Puy (ry )dry : 
0 R—rı+0(N) 


The second part integral has the estimation from A.2 


R R c*(t,d) g eat 
Dy vp (fv dr, = | 1 +o 1 e HNT , 1 +o 1))20 G 
E P(r) every PGE (1+0(1)) aa 


= o(l)e7o™, (3.2) 


The primary difficulty is to estimate the first part integral 


R—-ry,+@(N) 
[ Puyp(n)dry, (3.3) 


for which we need more fine analysis. We decompose it into two subparts. The first subpart has the 
estimation 


R-ry 7 = R-ry m P1 (Oxy) 
Puy (ry)dry = — ain r 1Ou yp (ry) dry (duy < Ty ry < R) 
0 o 40 tge "T 


R—ru T 
l f, ODP Our)dbuvp(r)dr 


(eT —e Tjd 


By -Nu 
= f O1 +o(1))20 4 am, 


= O(e776"), (3.4) 
For second subpart 
R-ry+@(N) R-ry+@(N) pr Oav 
| PurPlrs)dry = f LACON A (3.5) 
R-ry R-ry 0 fe —T 


we consider the angular integral p, at two cases since R — ru S ry S R—ru + o(N). 
(1) If ru < R—@(N), it has r, > @(N) . Together with the premier condition r, > @(N), we have the 
explicit distance expression from 2.1 


2 Buy Â, v 
duv =ru try + > msin ( 3 ) +0 (ér) , (3.6) 


hence the connection probability 


1 1 


duv R ` nNu+v-ZH . d ? 
l+e™T 1+Cye = sin? ($) 


for Buy <S Buy < 


m 
[ Pl (9) Puyd@ 


Buy 


T. So the angular integral 


N 
Puy 


II 


T ~ A 
pi(8)Puvd0+ f P1(9)puyd@ (Oy) > Buv) 


z T ai d-1 9 1 
= ofertu) f sin oo 7 dO 
o laai 14Cye ae sin? ($) 
n iyd-1 
(u+ -2 sin® 0 1 
= oferty (1of1) [E aea r io 
dl ae sint (>) 
95 T] sin? 1 
—(Nu+ -Z 
= o(e (nu+n H)) 4 1+o(1 Dl iks METER dé 
li {ene sint (5) 
95 sin’! @ T sin’! @ 1 
—_ —(Nu+v -Zp ) Po 1 f ) (| 
= o(e +o(1))O d0 )+(1+o0(1))0O R 
\+a+oaye( f° aol [aa — 
z o (eturman), 
utv-2 =l utny—2, 
where 05 is defined by sint (%8) = (e =) (Notice Cy = 1 +0(1) and e + ~~ > 1). The 
ô y 2 
last equality is from the following process. If @5 < 4, then 
B; cind-1 Ə; gd-1 
* sin O 40 =0( | a0) = 0 (6) = 0 ( eur -n) 
0 Ta o daa 5 
for (2%) < sin? (%8) < (%8), and if | < 05 < 7, then 
qs d-1 Os ind! m cind-l 
Z sin^ 0 6 sin®™ 0 sinf™ 0 
o< f ao< | d9, 
[ Ta 0 d,1 0 
so it is obvious Pa sedo —o(e emma), and there always has Ses in -e — a d0 = 
: eT sint (5) 


O| e~(m+m—4x) |) for0< O05 <1 


(2) If ru > R—a@(N), we luckily have the lower bound estimation of the distance from (1.4). Notice 

R-—3@(N) S ru — ry S duy S Fu + Ty, SO duy — œ as N — œ. We extend two sides of the hyperbolic 

distance formula (1.4) and get 

ebtut ory 
4 


efdur (1 4 e26dur) 
2 
so eb4uy = (1+ e 25duv))— LeorutSry (sin? Puy + 67 26% cos? Puy + 2S" cos? Puy + 2S lut) gin? ay, 
Hence 


((1—cos Oy) 4 et (] l cos Oy,y) + e725" (1 + cos Ou y) He 25") (1 — cos Bu v)), 


ef lduv-R) edur) -1 o$ lru trv—R) (sin? Sus pa uot us 4-28 cog? us q e 2 lrt) gin? Sus ) 


2 Puy +e 25" cos? 2 Sus 


(3.7) 


(1+ 0(1))e8+"-*) (sin? >? +67 25% cos? ioe eT? lut) sin 


ie 


Set 05 as before, i.e., sin? (8 )eS(utr-R) = 1. For 0 < 0,1 < 05 we have 


ef (duv—R) < 4(1 +o(1)), 


on the other hand, for 0 < 0,., < m we have 


duy—R rutty— A 0, 
eT > (1+0(1))e P sin 
Hence we have the upper bound estimation 
` 0 
Puy = | pit) d@ 
9 Ite -T 
65 0 m 0 
1 7 -e ao+ f 2 Ox dé 
0 Iper 9% 1+e`T 
95 sin’! @ T sind! @ 1 
[74+ ie n 
dl ô dl (1+o(1)je = sint(5) 


< e~ Mutt Bx) 
~ , 


and the lower bound estimation 


m m (6) fe] (6) in?! 0 7 
Buy = f P1(8) Puy dO = Pdo > | "Pit, do =0( | * S840) =o (eteina) 


l+e T 
(3.8) 


for the same process as before. 
Finally from the above two cases we generalize 


pd 
Pu,v = ‘i Pl (0) Puy d@ =0 Goa) . (3.9) 


Further for second subpart (3.5), we have 


R—-ry,+@(N) R—-ry+@(N) G 
i PuvP(Ty)dry = | O (e (wt 4H) 9 (r,)dr, 
R-ry R-ry 


R—rı+@(N) - 
ro) (| ettwa (7 dr, 
= G(e ™), (3.10) 


II 


Combining the estimation of two subparts (3.4)(3.5), we have the estimation for the first part integral 


(3.3) 
R—ry,+@(N) 
f PuyvP (ry)dry =0 (=) l 
0 


further together with the estimation of second part integral (3.2), we get 
7 R 
Prifuy = Uru, (JA = f PuvP (r,)dry =0 =) š 
0 


We may see that the connection probability Pr[f, = 1 |ru, 6, is mainly contributed by 0 < ry < R — ru + 
@(N) for ru > @(N). This property is different from the case of 20 > 1 in [16]. 


3.1.2) Privy = 1 ru, 6,,] on the condition ry < @(N). We discompose formula (3.1) into two parts as 
follows 


S R R—ru R 
Pr[Z v = Irv, Ou] = f Pu vp (ry )dry a i Pusplrv)dry + f PuyP (ty)dry. (3.11) 
R-ry 
For the first part integral, we have the estimation 


i oa P1(Oyv) 
0 


du,y—R 
09 14er 


R-ry 
| Puyp(r)dry dO, .P(rydry (duv <rutry SR) 
0 


R-ry oa 
Z n f, ODP Ouv)dduvplr dr 


Ru-Nu ( 29. Z — 26 D 
T O(1)(1 +0(1))20 4 — 5 an, 
= O(e 20M), 


For second part integral, since R—2@(N) < ry — ru S duv S Fu + fy, 80 du — œ as N — œ, we use the 
same analysis as second case in (3.5), there is 


Puy =0 (eum —Fn) ) 3 


so the integral 


R R m 0, , 
i Puvp(ry)dry = 1 PGi iird (3.12) 
R-ry R-r, J0 lte T- 
R 
= O(e (M+ Ft) y(r, dr, (3.13) 
R-ry 
= O(e7e™), (3.14) 


Together with two parts, we have 


D R 
Prila = Uru, 8u] = f Puvp (rv)dry = O(e 22), 
0 


3.1.3 Privy = L|ru, 6, for any ry. From §3.1.1 and §3.1.2, we have the first conclusion (1.10) in 
Theorem 1.1, 


: R 
Priluy = Uru, 6,,] = i Pu,vP (r,)dry = O (e?u) 
f 0 


for any r,. Figure 1 shows the simulation results which is consistent with our analysis. On the other 
hand, we also get a fact: for O < ry < R — ru, Puy = O(1); for R — ru S ry SR, Puy = O (e~ (Mutu), 


3.2 Mean degree 


By the simple computing, we have second conclusion (1.11) in Theorem 1.1, 


(k) = (N— 1) ['pjete2™) ar =0 ((v- Hye?" ) 
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chinaXiv 


log(Pr) 


FIG. 1. Plot comes from simulation of Pr[J. = 1|ru, 6,] by the law of large numbers. The vertical axis is log(Pr) and the abscissa 
axis is radius r,. The data points with N = 2!? from lower to upper take parameters d = 3,20 = 1/2,t = 1/2; d= 1,20 
2/3,t = 1/3; d = 1,20 = 1/2,t = 1/2 respectively, and € = 1. Comparison of three almost straight lines is consistent with 
equation (1.10) 


3.3 The probability of k-degree 


Next we analyze the probability of D, = k, where k is a positive integer. Let constant A € (0,1). Then 


Pr[D, =k] 
a ') (Prila = 1ra) (Pella = 1ra) p Cru)dru 
pn NT) Erlar = trl) Petes = Ural)" peara 
= f E ') (Prusy = 1rul)* (1 — Prusy = 1ra] E p (ru)dra+0(1). 


Now to analyze the first integral, if let £ = Pr{/,,., = 1|r,] = © (e720) and we can see t + 0 uniformly 
for ÀR <r, < R, and 


dai _ (1-+0(1))expl-(W — 1) x expl- zi Eya 


0<€ <t, o(1) = O(e 24) when ru > AR. So 
LE ') (ri (Prlfigy = Urel)€ (1 —Prlfer = Ural)¥ E pedr 
_ LC : ra Pi ini oi ap Wai aia 
exp [sar (Plan = Ura)" odra 
1 


R k 
EE I p gg (O7 D Prka = Ural) exp [-(N = 1) (Prlluy = Hea) 


N-1 
x exp -a (Priliy = id) P(ry)dru- 


Since 
(N —1)@(e72°%#) < (N—1)PrlIyy = 1|ru, 61] = (N — 1)O(e720™) < (N— 1) 0 (e204 #41), 
xk 7 


when take 2y = InX, then (N — 1) Pr[uy = 1|ru, 6] — œ as N —> œ, and since = 
for x > k, then 


E ((v- 1)Pr[luv = 1 rus 6,]) exp |-w- 1) (Pilles = Ire 8d) +0. 


So we get 


R /N—] a4 Nk _, \N-1-k 
| ( i ) (Prt = tr il) (1 = Prifuy = Uru, 4.1) 0(ru)dry > 0, 
AR É 


as N — oo, Finally we have third conclusion in Theorem 1.1, 
Pr[D, =k] > 0 


as N — œ for any given k > 0. Simulation (figure 2) also shows the same tendency. 


4. Expectation of incomplete and complete triangles 


In this section we will compute longer link and prove Theorem 1.2. 


4.1 Incomplete triangle for the arbitrary dimension 


In this subsection, we compute the probability P(A(u,v;w)) of the event A(u,v;w). First we define 
the angular integral Pu v;w by a conditional probability which describe that nodes u,v at the fixed radial 
coordinate r„,ry connect to a node w at the fixed radius rą. We pivoted the radial direction of vertex 
w as the center axis and the connections of {u,w} and {v,w} are independent events, we may compute 
Pu.v;w by integrating over the relative angles 0, w, 0, of u,v about w respectively(figure 3). 


T pt 
Puy;w = Pu,wPv,w = f f pı (Ou w)P1 (Oy w)Pu,wPvw dOu wd Oy. (4.1) 
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FIG. 2. Plot of probabilities of node degrees having various values of k, against the network size N, for d = 1, T = 1/2, 20 = 1/2, 
v = l and € = 1. one can see how Pr(k) — 0 for fixed k as N > œ. 


FIG. 3. Relative angle of the triangle in d = 1 and d = 2 


Notice for d = 1, we have py w = Io Paw dOu w since the angular symmetry for the hyperbolic distance 
duw. SO we have 


R R R 
Pauw) = f S h Berpe) ldradr dr, 


R pR pR 
f f f birbaP edP (rw)drudrsdr 


R 
z f PPa = ij Prie = 1rwldrw 
f | 


R 
[ © (h) P(rw)drw 
0 


Ry (eTe Tad 
—40nw 
1 o (e on a +0(1))26< 5" dy 


(0) Gas . 


Hence we get the formula (1.12) in Theorem 1.2, 


E(A)= (3ra ames 3 )e =). 


4.2 Complete triangle for the arbitrary dimension 


II 


In this subsection, we consider complete triangles and finish the proof of second part in Theorem 1.2. 
If we choose triangle formed by (u,v,w) pivoted at w, for d = 1, there is the similar angular integral 
definition 


P T pt 11 
Puy,w = P(A(u, v,w) |ru, Fy, rw) = f i PuwPvwPuy 7d Pu wd w» (4.2) 


where Ou = |Ou,w — Ow] or min{ Oy y+ O),27 — (Ou,w + A) } (refer ro A in figure 3); for d > 1, the 
angular integral 


m T T 
Pu,v,w = P(A(u, V, w) = Uru, Fy, rw) = f f f Pu,wPv,wPuyvP1 (Ow) Pt (Ow) P2(O,2)dOp wdO) dO, 2 


where 


pi(01) = sin?! 0) Ian, p2(02) = sint? 02 /Ta2 


and 0,2 is the angle between uw-plane and vw-plane (refer to B in figure 3). It is obviously pyyw < 
PuwPvw- The probability of the event A(u, v,w) is 


R rR pR 
Pw »W)) = f f h Pune (rade OP )drdrdry. 
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If rw > @(N), there is 


P(A(u,v,w)|rw > @(N)) 


In the following, we consider P(A(u,v,w)|ry < @(N)). We start from the range rw < @(N), fu, ry S 


R R pR 
f fff Porp edP re (rw drudrvdry 
N) Jo Jo 


IN 


R R pR 
f f f Pu,wPv,wp (Tu)? (rv) P (rw)drudrydrw 
N) Jo Jo 


R R R 
= 1 prw) f PuwP(ruddra | P (ry) Py wdrydrw 
N) 0 0 


R 


= (0) Gua P(rw)dry 
N) 


B 
/ H etongge2(m-4l dy) 
o (N) 


= o(P(A(u,v;w))). (4. 


& 


II 
© 


and there is Puw = O(1) since Puw, Prw, Pu,v are O(1), so 


on the other hand, we have 


IN 


a 


P(A(u,v,w)|Tw < O(N), Tasty S z) (4.4) 


oN) p pł 
S S S Pune rade de (rw )drudrvdry 
0 0 0 
oN) pf p3 
f f I Olip@ pt, )pt drar 
o(N 
off" m ar f p(ry )dr,p( ry) )dry, 
0 
o(N 
og ns 
0 
o (etone 200 (N ) 


R 
P(A(u,v,w)|rw S o(N), 5 Sru SR, rn <x) (4.5) 


NI 


a(N) fR p8 
1 ii PuywP (Tu) P (tv) P (tw)drudrydry 
2 


oN) fR p8 
f n f P(ry)P (Tu) P (rw)drudrydrw 
2, 
@(N) fR R 
f J peddi [ alr) dryp (ry dis 
0 k 0 


@(N) 
f e #4 p (ry)drw 


0 
—302y 20) (N) 
e ? 


e 


and 


R R g 
P(A (u,v,w)|rw < O(N), D <r, SR, ru S 5) < gerne). (4.6) 


and for the triangle at the below situation, we choose to the triangle pivoting at v 


R 
P(A(u,v,w)|Tw < @(N), 5 


@(N R pR 
= f J hi PuyvwP (ru)P (ry)P (rw)drudrydrw 
0 743 


S ru,ry SR) (4.7) 


a(N R pR 
= f J j: PuwyP (Fu)P (rv)P (rw )drudrydrw 
0 RJR 
pivoting at v @(N) PR PR 
i [ I Puy PwoP (Tu) P (Tv) P (rw )drydr,dry 
Z“ z 


R R 
| (J, Pu vP (ddr) pP(ry)dryp(rw)drw 
2 


/N 
z 
& 
2 


R 
2 


R 
e270 5 (ry )dryp (Tw )dry 


LA 
o— 
e 
= 
ea 


ai o —20fH 
< Bye P(ry)dry 
0 
&, 


< eo FOF p200 (N). 


H 
Together (4.4),(4.5),(4.6), and (4.7), we have 
P(A (u,v, w)|rw < @(N)) = o (P(A (u,v;w))). (4.8) 
Hence from (4.3) and (4.8), we have the order estimation (1.13) in Theorem 1.2, 
P(A(u,v,w)) = o(P(A (u,v;w))),E(T) = 0(E(4)). 


Next, we will give the lower bound estimation for P(A (u, v,w)) > e7398 , which will deduce the lower 
bound estimation (1.13) in Theorem 1.2, 


N 4 
E(T)= (; PAu w))2 ( jee (4.9) 
From the similar process of (4.4), we get 


P(Atu, V, w)|Fw, Fufy < 


) 


R ph fh 
2 [ [ [ DuvwP Fa) P (JP (ie ldrudrdry 


- 1 2 f 2 f O T E A 


(2) eo), (4.10) 


on the other hand, we consider 


R R 
P(A(u,y, wlz SWL Rn aS z 
Rod f3 
— I. f f Pu,v,wP (ru) P (rv)p (rw)drudrydrw 
2: 


Rot ph 
= Í i I Pu wPvwO(1)p(ry)P (Tu) P (rw )drudrydrw 
2, 
R R R 
(0) (/ Į Pu wP ry Į Pv.wpP ejan Pirin) 
2 
ù Aoi 3 R-ry 8 
= 0 | f a PuwP (Tu)dru [ F PyvwP(ry)dry| P (rw)drw 
R 0 R—rw 0 R-n, 


R R—rw R : 
= l Į O(1)p(ru)dru + ? oleonapo )ar 
0 


R-ry 


II 


R—rw R 
| O(1)p(r dr, + |? ole Cna p(n) Ptr) (4.11) 
0 


R-ry 


R 
= o( A Oe 2°™)9(e2°™)p(r)drv) 
2 


= O(c 30%), (4.12) 


where third equality from last (4.11) involved the below details: 
R-ry Ru-—Nw (ea "u e plu 7 Z 
[ i a (/ 2202y dn, | = 0 (e °°™ f(nw)) 


where f(y) = e720(2u—Mw) [Ze (en = ed 1)4dn, and 


2 - (wtu -Zy) p “(tne -&y) (Co ea md -20 
A O (e Nwt+Nu-AH )p(ru)dry =@0 m e Nw+HNu— AH Zoia dNu =0 (e ™ 9(nw)) ; 
TEW AH '!lw 


G Zy 20 20 
where g(nw) = e(1-20)(#u—Nw) Za- m eM (e Tu —e~ “a ™)4dn,,, and by the general monotonic method, 
we can get f (Nw) +g8(Nw) is initially increasing then decreasing when Nw varies from Zu to Zy, so we 
have f(nw) +8(w) = min{ f(#“),¢(#u)} > co(d,o), where the constant co(d,o) is only depending 
on c,d for enough large N, on the other hand, the upper bound of f(1,,) +g8(nw) is obvious, so we have 


R-ry R 
j O(1)p(ru)dry+ |? O(e MBit) (ry dry = O (e?) 
0 R-ry 


Combining (4.10) and (4.12), we have 


R ph of : 
S S f Pan rupov)plrw)drudrrdry =0 (23%) (4.13) 
0 0 0 


hence P(A (u,v, w)) Z e~3°##, Question: P(A(u,v,w)) = © (e7924)? 


FIG. 4. Two Complete Triangles 


5. Variance estimation of incomplete and complete triangles 


We first estimate E (T?) and E(A”). We need to set out the specific expression Zy = In ¥. 


5.1 Estimation of E(T*) 
LEMMA 5.1 E(T?) = (1 +0(1))E?(T). 


We mark two complete triangle events A(u1,v1;w1), A(u2, v2; w2) formed by two groups of vertexes 
U1,V1,W1 and uz, V2, W2 respectively. First we have 


2 
A tag oe 3 ( y Lato =E(T)+ Elta nan Atos): 6-1) 
(u,v,w) (u1,v1;w1)#(u2,v2;w2) 
for second part events, there are the below three cases (figure 4). 
1. {u1, v1, w1} O {u2, v2, w2} = 0. 
2. {u1, v1, w1}, {u2,v2,w2} have only one common vertex. 
3. {u1,v1, wi}, {u2,v2,w2} have exactly two common vertexes, i.e., they share a same edge. 


We denote the event by T; corresponding i-th case, i = 1,2,3. 
For Case | (refer to A in figure 4), there are (3) possible choice for {u1,v1,w1} and leave CS) possible 
choice for {u2,v2,w2}. So 


E(Ti) = 3 (1 PAv wP Alw) = (1+o(1))E7(T). 


For Case 2 (refer to B in figure 4), there are et ) possible choice for {u1,v1,w1} and if w; is common 
vertex, there leave io 3) possible choice for {u2,v2}. 


P(A(u1,vi,w1), A(u2,v2,w1)) 


R fR fR fR fR 
0 f f f Pun vi wi Puz,va wi P (Fu, )P (Tvi )P (Fu )P (Tva )P (Fw; dru, drv dru dry drw, 


IN 


j 

R fR fR fR fR 

LSS | baron Beroe Bi anPi 2 PEP n)a) (Pag) dr dy dry Arndt, 
OB e 0%, 


/N 


so considering the lower bound of E(T) in (4.9), we have 


E(Tp) 2 


II 
(09) 
P 
w z 
No Ne 


< 3 @ @ = P) oaie to 
3/\ 2 
2 
= o L ae 


= o(E*(T)). 


(5°) Ptastan ven) Azv) 


For Case 3 (refer to C in figure 4), there are Cy possible choice for {u;,v1,w1}, if v,,w, are common 


vertexes, there leave es O possible choice for up. 


P(A(u1,¥1,W1), A(u2,V1,W1)) 


< 
< Oke tn, 
E(T3) = (ZC T PAv), Avm) 
(UT neen 
= o(E*(T)) 


Combining all cases, we get lemma 5.1. 


5.2 Estimation of E(A”) 
LEMMA 5.2 E(A7) = (1 +0(1))E? (A). 


Similar, we have 


2 
E(A*)=E ( y Lair) =E(A)+ 


(u,v:w) (u1 ,v1;w1)#(u2,v2;w2) 


R fR fR fR 
LL Lf Panbrnbump ann) En)dradradradrna 


E (14A (ui vi iwi )} HA (2,232) })- (5.2) 


We may classify second part events into the following eight cases (figure 5). 


1. {u1,v1,wi}M {u2,v2,w2} =0. 


The following case is that {u1, v1, w1 }, {u2,v2,w2} have only one common vertex: 


2. wy =W2, {41,1} O {u,v} = 0. 


3. uy =u2, {v1, wi} A{v2,w2} = 0. 
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FIG. 5. Two Incomplete Triangles 


4. u = w2,{v1,wi } N {u2, v2} = 0. 
The below case that {u1, v1, w1 }, {u2,v2,w2} have exactly two common vertexes: 


5: uy u2,V1 v2,W] w2. 


6. ui = w2,V1 = V2,W1 Æ u2. 
7. Uy = U2,W] =W2,V1 É v2. 


8. uy = w2,w1 = U2, V1 E v2. 


We denote the event by A; corresponding i-th case, i = 1,2,--- ,8. 
For case 1 (refer to I in figure 5), there are 3(7) possible choice for {u1,v1,w1} and leave a0 5 >) 
possible choice for {u,v2,w2}, so 


BAS (3) ("5 a Pie Gawd oA). 


For case 2 (refer to II in figure 5), there are 3 E ) possible choice for {u1 , vı, w1 } and leave (Y z) possible 


choice for {u2,v2}, so 
P(A (u, v1;w1)A (u2, v2; W1)) 
R fR fR fR fR 
= SS Bambam Brow Broo PEP a)n )dradrndradrdr 


R R R R R 
[ i Pum PCa )drn, f Prion P (dr, f Pemb Can) | Pram P addr} Blt dro 


R 
O (i ptr )dry, 


Ay (ea -e IF wi \e 
= —80Nw 
z f o € ) (1 +0(1))20< — dnw, 


= @ Gus 


and 


E(Ay) = 3(3) (“;*e ~ = o(E?(A)). 


S For case 3 (refer to II in figure 5), there are 3 (3%) possible choice for {u1,v1,w1} and leave Ce) 
possible choice for {v2,w2}, so there are totally 6(4 ) Cy 3 ) choices for such case. 


Pe vi;w1)A (u1, v2;Ww2)) 


[i T [ f Êv, wi Pwy ay Puy wa Pwa, vaP (Fu, )P (fv; JP (Fwi P (fva )P (Fws Jdr u, dry drw, dry, drw, 


R 
= f ii Py, w in Puy wi ie Puy wo in p(ra)Pramdr| pre )dra| pra arn pm dr | p (ry, jaty 


II 


R R 
ie w if Puy wi if Pum O (ema) plrn)drns P (Tu Jaru | P (Tw, jarn | P (rv )dry, 


R 
aT w if Buy w (Bue °° p (ru, Narn P (rw; Jar | P (rv )dry, 


IN 


R : 
rs "an w OB e’? O (e701 )p (fwi jarn | P (rv )dry 
0 


R 


IN 


ne 2°44 O (By)? p (ry, dry 


where we use 7 
f Puy Gay p(rv)dry < O(Ru)e 7°". 
0 
So 


BAA = 6(3 ) @ E °’) O(R2 jea = o(E?(A)). 


For case 4 (refer to IV in figure 5), there are 3(7) possible choice for {u1,v1,wi} and leave Ce) 


possible choice for {u2,v2}, so there are totally 6() (3°) choices for such case. 
P(A (u1, v1; w1)A (u2,v231)) 
R pR pR pR pR 
= f f [ i f Pu wi Pvy wi Buo uy Pwa, P Tay P (Fv P (Fwi )P (Fv )P (Fua dru dry drw, dry d'u 


R R R R R 
f if Purn P(n ir | if Puy i Pn Plinth] i Pera P (Pra) | PUP brn | Plt) 


R R 
= | O (e7? ) / Pum O(n J (6-2 )p (Fn Jr P (Tw, )drw, 
0 0 


IN 


[9 (e291 op (Fy) 
< O(Bye +O") 
where we use k 
[Burm @(e20% JO(e-2° pra, dra, < O(€ 2%), 
So we have 
E(Ag) =6 n ) (* A `) O( Bye 404") = o (E2(A)). 
For case 5 (refer to V in figure 5), there are 3%) possible choice for {u1,vı,w1} and leave O) 
possible choice for w2, so there are totally 3 ie ) ( ae ) choices for such case. 
P(A (u1, vi;w1)A (u1, v1;w2)) 


ES L [i ers Bean Page, 0a) )0 0) Fg dre dra ry 


R 
[ H o Piwi ia Buy pP (teres | Pra dria | Bm re | BCP 


R 


R 

< fle Pn One®") | 0 
R 

< fÉ ARa ocem )p lr, jdr 
0 

< (Bhen), 


So we have 
E(As) = 3(3) (T Pata vise). vwa) < 3(3) Ci `) O(#ye 494") = o (E?(A)). 


For case 6 (refer to VI in figure 5), there are aC) possible choice for {u1,v1,wi} and leave Ca) 
possible choice for w2, so there are totally 3 (3) ( oa) choices for such case. 


R R 
PEACH, vi w)A ua vm)) < f f f f Pvy wi Pwi ey Pusu (Tuy )P (Tr, )P (Fwi )P (Fus ) Arup dru drw, dry, < O(Rye 9n), 


So we have 


E(A6) = 3(%) ra °’) P(A (u, v1;w1)A (u2,v1;u1)) < 3(%) Ca) O( Rye 40%#) = 0 (E*(A)). 


For case 7 (refer to VII in figure 5), there are aC ) possible choice for {u1,v1,w1} and leave w ‘i 


possible choice for v2, so there are totally 3 (3 ) Oe ) choices for such case. 


P(A(u1,V15;W1)A(u1,¥25W1)) 


< SS LF ina Pus Prom P Fu )P Fo )P Pe) P lw, Jr ro dren 
- f | : “Bamba dr | | I "Pray Pla dra | [ “Pram Prada} Plt drw 
= fi om 0 (2 JO(e2™ Jp (Fw, dr 
= O(c 204), 
So we have 
E(Aq) =3(3) (87?) Paw nna) < 3(5) ow O(c 20H) = o (E2(A)). 


For case 8 (refer to VIII in figure 5), there are 3%) possible choice for {u1, v1, w1 } and leave C 
possible choice for v2, so there are totally 6X ) Ce à ) choices for such case. 
A (u1, v1;w1)A (w1,v2;u1)) 
T [ [ a Pv wi Pwi u Puy vaP (Tui )P (Tvi )P (Fwi )P (fva dry, dru drw, dry, 


< O (Bye Aoir), 


So we have 
ra= (y) (7?) Aniwa <6(3) ("> ?) 02k oa). 


Together with all cases, we have lemma 5.2. 


5.3 Variance and Chebyshev’s inequality 


From previous subsections, we have 
D*(A) = E(|A —E(A)|*) = E(A*) —E°(A) = o(E’(A)), 


and 
D*(T) = E(|T —E(T)|) = E(T°) -E° (T) = o(E?°(T)), 


if Zy = ln y, Further by Chebyshev’s inequality 


P(IX - E(X)| > £) < 


if we have £E? (X) in steady of £, then we get 


(PREA e)s 


d=1,20=1/2,1=1/2,v=1 


§ 4 (2*10,0.1868) 
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global clustering coefficient 


FIG. 6. Plot of tendency of global clustering coefficient against the network size N, ford = 1,20 = 1/2, t = 1/2, v = l and ¢ = 1. 
One can see the tendency of C2 (G) > 0 as N > œ. 


e., A =E(A)(1+0,(1)). 


Similar, 
T =E(T)(1+0,(1)). 
6. Global clustering coefficient 


The global clustering coefficient C2(G) of a graph G is defined as 


So from previous section, we have Theorem 1.3 


_3T(G) _ 3E(T)(1+op(1)) p 


CG) = @ Aro] 


+O, 
as N — œ. Simulation result (refer to figure 6 ) also supports this conclusion. 


7. Conclusions 


In this article, we give a precise characterization of degree and clustering in the random geometric graph 
on the hyperbolic ball of any dimension for the small parameter t < 1,20 < 1 and the given relation 
of By = In%. We get the to-zero tendency for the k-degree and the global clustering coefficient as 
N — œ. Maybe it doesn’t meet the desire for the power law of small parameter in the real world, but the 
interaction between the density parameter o and the region expanding speed is reflected by this article 
and the previous article [16], which reveal that mobility may has an important significance as a generic 
mechanism in networks. Some techniques and analysis method are firstly used in this article and we 
hope that our results and methods may be applied in machine learning since hyperbolic geometry has 


been widely used in this field for the quasi-isometric embedding proof of Sarkar [8, 15]. Here we use the 
global clustering coefficient to describe the clustering phenomenon, and there also other definitions for 
the clustering tendency, so it is interesting to try other clustering definition for this model. We recently 
know that Fountoulakis give the finer characterization about clustering in [13]. It is a natural step to 
follow their work for the case T < 1,20 < 1. There are many other models of hyperbolic geometry, for 


example in [6, 7], and some researcher advices us that Klein model in [6] is 


well-suited for calculation, 


so it is a good choice to apply our method to such model. We also think that the random geometric graph 
model on the hyperbolic ball can analyze the social network since there exists some common principle 


as we mentioned in Introduction. 
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A. 
For ru +ry —R > @(N) we divide the integral (1.5) as follows 
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and the last equation is from the estimation J = 6 fọ sin’ “8 48 — ¢*(¢,d)5* + 0(5*) in the Appendix 


ô+sinT ($) 
Z 2y—-Nu-MWw 
A.1 of [16] with ô =e = ~. Finally, we have 
5 4 c*(t,d) Rr —Nu-— Nv 
Buy =} P1(8)Puyd = (1401) eo, (A.2) 
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